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1  Introduction 


Research  conducted  o\  er  the  past  decade 
has  led  to  the  dev  elopment  of  a  number  of  innov  a- 
tive  procedures  for  estimating  the  >  ields  of  under¬ 
ground  nuclear  explosions  through  s>stematic 
analyses  of  digital  seismic  data  recorded  from 
these  tests.  In  addition,  a  wide  varietj  of  new  data 
regarding  the  geophysical  env  ironments  at  Sov iet 
test  locations  have  now  become  available  as  a 
result  of  the  Joint  Verification  Experiment  (JVE) 
and  associated  data  exchanges.  The  objective  of 
the  research  program  described  in  this  report  is  to 
integrate  all  these  new  capabilities  and  data  into  a 
comprehensive,  prototype  system  which  can  be 
used  to  derive  optimum  seismic  estimates  of  ex¬ 
plosion  y  ield.  More  specifically ,  it  is  to  implement 
a  fiexible  interactive  software  system  in  which 
y  ield  estimates  based  on  a  w  ide  v  ariety  of  different 
seismic  magnitude  measurements  can  be  eff  ^ciently 
determined,  merged  vv  ith  all  av  ailablc  information 
regarding  the  test  location  under  consideration  and 
statistically  combined  to  obtain  a  unified  seismic 
estim.ate  of  explosion  y  ield  and  quantitative  mea¬ 
sures  of  the  uncertainty  in  that  estimate. 

A  preliminary  prototype  version  of  a  sys¬ 
tem  designed  to  achieve  the  above  objectives, 
designated  the  Yield  Estimation  System  (YES), 
has  been  implemented  in  a  Sun  color  workstation 
(SPARES tation)  env  ironment  at  the  DARPA  Cen¬ 
ter  for  Seismic  Studies  (CSS)  using  software  built 
upon  the  framework  of  the  X  Window  graphics 
and  Oracle  database  management  .systems.  This 
initial  version  focuses  on  explosions  at  the  Sov  ict 


Shagan  River  test  site  and  on  y  ield  estimates  ba.sed 
on  five  different  seismic  magnitude  measures  de¬ 
rived  from  teleseismic  P  and  surface  wave  tmb, 
Mpspei.,  Mo)  and  regional  L»  and  P  wave  (Ml*;. 
Mpn)  data.  Available  test  site  information  is  pre- 

T\i 

sented  to  the  analyst  in  the  context  of  a  SPOT 
satellite  image  of  the  Shagan  River  region,  in  a 
format  which  permits  the  analyst  to  interact  digi¬ 
tally  with  the  image  to  easily  extract  and  display 
information  regarding  the  explosion  .source  envi¬ 
ronment.  A  waveform  database  consisting  of  more 
than  10000  digital  seismograms  recorded  from 
Shagan  River  explosions  at  stations  of  the  GDSN, 
USAEDS,  NORSAR,  CDSN  and  IRIS  networks 
has  been  assembled  for  this  sy  stem.  The  graphical 
user  interface  to  this  system  is  completely  menu- 
driven  and  moiise-activ  ated  and  has  been  designed 
so  that  no  keyboard  entry  is  required  of  the  opera¬ 
tor.  Thus,  the  seismic  data  can  be  directly  input  to 
the  V  arious  datu processing  modules  using  asin;ple 
menu  .selection  procedure  to  obtain  the  magnitude 
measures  and  associated  y  ield  estimates. 

This  report  presents  a  summary  of  the 
current  status  of  the  ongoing  research  investiga¬ 
tions  directed  toward  the  development  of  an  im¬ 
proved  seismic  yield  estimation  capability  for  un¬ 
derground  explosions.  The  system  design  criteria 
are  reviewed  in  Section  2,  where  the  characteris¬ 
tics  of  the  data,  analysis  tools,  databa.se  relations 
and  graphical  usci  interface  arc  described  in  the 
context  of  their  integration  into  a  comprehensive 
software  system.  This  is  followed  in  Section  3  by 


TM  SPOT  data  arc  copyriglucd  by  CNES  { 1986.1987). 
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an  overview  of  system  capabilities  in  which  the 
functionalit)  of  the  current  protot>pe  sjstem  is 
graphicall}  illustrated  using  display  s  of  the  s«-reens 
encountered  by  an  analyst  in  a  typical  processing 


session  for  aselectedij.e.,  the  JVE;  explosion.  The 
report  concludes  with  Section  4  which  contains  a 
sammarj  and  a  discussion  of  future  plans  for  the 
development  of  the  YES. 


2  YES  Software  Design 


The  '^'ES  is  designed  to  determine  an  opti¬ 
mum  explosion  yield,  the  uncertainty  in  this  yield, 
and  a  statistical  assessmeiil  of  the  results,  using  a 
wide  variety  o*"  analysis  tools  applied  to  .seismic 
signals  generated  by  under  ground  nuclear  explo¬ 
sions.  Such  a  s>  siem  requires  .i  flexible  design  that 
allows  analysis  tools  to  be  easily  added  and  yet 
allow  s  all  or  part  of  die  analy  sis  tools  and  al  1  or  part 
of  the  data  to  be  used  for  my  gi\  en  y  ield  estimate. 
Furthermore,  bccai'se  the  yield  estimates  depend 
on  seseral  types  of  analysis  and  a  variety  of  data 
which  will  not  only  vary  from  one  event  to  the 
next,  but  vv  ill  also  change  as  new  data  becomes 
av  ailable  for  a  single  event,  the  system  must  allow 
identification  of  all  of  the  data  and  proce.ssing 
steps  that  correspond  to  any  yield  estimate  made 
with  the  sy.stem. 

In  this  section  of  this  report,  we  di.scuss  the 
YES  softw  are  design,  the  conceptual  model  for  the 
system,  the  requirements  and  constraints  of  the 
system,  how  the  database  has  been  designed  to 
handle  general  analysis  tools  and  how  the  system 
has  been  structured  to  achieve  design  goals.  This 
is  intended  to  prov  ide  a  description  of  the  undcrly 
ing  structure  of  the  system  that  forms  a  platform 
for  operation  of  YES  modules. 

2 . 1  Conceptual  Model 

YES  is  an  example  of  a  general  class  of 
.systems  built  on  four  basic  elements. 


®  data 

•  analysis  tools 

•  a  database  management  sy  stem 

•  a  graphical  user  interface. 

In  this  case  the  data  include  .seismic  wave¬ 
form  data  from  underground  nuclear  explosions, 
related  data  such  as  event  locations  and  material 
properties  at  test  sites,  parametric  data  such  as 
magnitudes  derived  from  seismic  data,  historical 
data,  and  image  data.  Tire  analy  sis  tools  operate  on 
these  data  and  extract  and/or  display  information. 
Examples  of  analy  sis  tools  arc  programs  to  display 
waveforms,  measure  magnitudes,  and  ov  erlay  event 
and  geographic  infonnation  on  maps  and  images. 
Tiic  data  and  all  of  the  analysis  results  up  to  and 
including  the  inferred  explo.sion  yield  are  main¬ 
tained  within  a  relational  database  management 
.system.  This  system  performs  both  the  function  of 
keeping  track  of  raw  data  and  related  information 
and  the  function  of  maintaining  a  procc,ssing 
history  of  each  y  ield  estimate.  The  graphical  user 
interface  allows  the  user  to  interact  with  the  data, 
database,  and  analysis  tools.  All  user  interaction 
takes  place  through  the  graphical  user  interface. 
For  YES.  this  interface  has  been  designed  to  mini¬ 
mize  the  amount  of  keyboard  entry'  and  to  allow 
nearly  all  .selection  to  occur  through  mou.se  driven 
menu.s  and  display.s. 

YES  consists  of  a  hierarchy  of  programs. 
.At  the  top  level  i.s  a  master  program  vvho.se  primary' 
function  is  to  start  the  other  analysis  programs. 
The  lop  level  program  allows  the  user  to  select  the 


test  site,  e\ent,  and  phase  to  be  processed,  and  the 
anal)  sis  tool  to  appi)  to  the  data  corresponding  to 
that  selection.  It  then  starts  the  anahsis  program 
with  the  proper  input.  The  master  program  starts 
some  processes  automaticjll) ,  experience  has 
shown  that  user  time  can  be  saved  bv  running 
certain  coniputationall)  -intensive  modules  as  .soon 
as  input  data  arc  read)  for  them.  B\  design,  the 
master  program  is  both  small  and  simple.  Because 
of  this,  it  is  easv  to  expand  the  program  to  add 
additional  anahsis  modules  as  required. 

A  difficult  problem  in  the  design  of  the 
inter  _tive  s).stem  was  to  make  it  robust  and  ea.s) 
to  use  for  an  operator  w  ith  limited  training  while 
retaining  vcrsatilit)  .sufficient  for  an  expert  to 
appl)  the  .s)  .stem  to  complex  c.Lses.  Our  .solution  to 
this  problem  h.is  been  to  provide  user  control  over 
operational  characteristics  of  the  anal)  sis  modules 
to  allow  performance  of  complex  operations  but 
also  to  build  in  default  settings  that  arc  appropriate 
under  most  circumstances.  Thus,  a  user  can  oper¬ 
ate  the  s\  stem  in  a  routine  w  a  v  bv  simply  follow  ing 
a  well-defined  set  of  steps  from  addition  of  new 
data  through  estimation  of  explosion  )  icid.  How 
ever,  if  special  problems  occur,  such  .is  an  anoma 
lous  result  or  the  need  to  modify  a  )  ield  estinnile  in 
response  to  other  know  ledge,  an  expert  can  adju.si 
the  operation  of  the  modules.  In  either  case,  the 
chosen  parameters  arc  stored  in  the  datahese  .ind 
can  be  recovered  together  with  the  resulting  yield 
estimate  at  any  time. 

2,2  Implementation  Framework 


The  graphical  user  interface  of  the  Yield 
Estimation  S)stem  is  b.i.scd  on  the  X  Window 


Sv.stem.  X  was  designed  specificaii)  to  allow 
hardware  independence,  to  foster  cmsc  in  porting 
applications  to  machines  other  than  tho.se  for  w  hich 
they  were  developed,  and  to  permit  running  of 
applications  on  one  computer  while  displaying 
their  output  on  another,  even  if  the  computers  arc 
of  different  manufacture.  YES  was  designed  and 
implemented  on  Sun  SPARC  computers  le.g.. 
SPARCStationst.  but  transfer  to  an)  other  .system 
that  supports  X  and  UNIX  could  be  accumpILshcd 
easily  .  The  X  interface  luis  beer,  w  ritten  using  the 
X  toolkit.  The  X  toolkit  enforces  an  object-ori¬ 
ented  approach  to  programming  by  ..ombining  the 
windows  and  the  operations  on  the  windows  into 
"widgcLs".  YES  uses  .several  widget  sets  including 
the  Motif  widget  .set  from  the  Open  Software 
Foundation,  the  X  widget  .set  from  Hewlett  Packard, 
the  Athena  w  idget  .set  from  MIT.  graphics  widgets 
written  by  Telcdyne  Gcotcch  and  SAIC.  and  spe¬ 
cial  purpo.se  w  idgets  for  YES  written  by  S  -CUBED. 

YES  is  written  in  the  C  and  FORTRAN 
programming  languages.  C  was  chosen  bccau.se 
programmers'  calls  to  the  X  Window  System  pro 
cedures  are  in  C.  and  becaaseC  is  w  ell  .suited  to  t..c 
design  of  complex  .systems  with  a  variety  of  data 
.structures.  FORTR.^N  was  used  bccaasc  it  is 
optimal  for  certain  types  of  <.ompuIatiunal  anal)  - 
sis.  and  necau.se  iLs  use  permitted  the  inclusion  of 
many  previously  existing  analysis  tools  into  YES. 
As  a  result,  the  top  level  programs,  interactive 
modules,  database  interface,  and  graphics  routines 
arc  written  in  C.  and  when  appropriate  these 
routines  call  FORTR.-XN  subroutines, 

YES  U.SCS  the  Oracle  relational  database 
management  system  as  its  database  manager. 
Oracle  was  cho.sen  because  it  is  u.scd  as  the  data¬ 
base  manager  at  the  Center  for  Seismic  Studies. 


and  because  it  is  available  for  a  wiJi,  variety  of 
hardware  platforms.  YES  uses  a  library  of  proce¬ 
dures  to  access  the  database,  and  all  Oracle-spe¬ 
cific  calls  have  been  isolated  m  that  they  can  be 
replaced  with  a  minimum  of  effort  if  necessary. 
The  relational  database  tables  are  independent  of 
the  particular  database  management  s>stem  used 
to  store  and  access  the  tables.  YES  can  be  config¬ 
ured  to  use  an  internal  database  manager  if  Oracle 
is  not  available. 

All  database  tables  are  documented  in  a 
standard  format  before  being  created  in  Oracle  or 
used  within  the  system.  All  Oracle  creation  and 
load  scripts  as  well  as  the  header  files  that  define 
internal  data  structures  and  external  file  formats 
are  generated  directly  from  the  documentation. 
This  way  t>pographic  errors  are  avoided  and  any 
changes  in  the  database  structure  can  be  made 
throughout  the  system  by  changing  only  the  docu¬ 
mentation  and  regenerating  dependent  files. 

2 . 3  Database  Structure 


The  Yield  Estimation  System  is  constructed 
on  a  substrate  of  a  database  management  system 
and  a  database,  which  serve  to  manage  all  data 
required  for  communication  between  modules. 
Each  module  extracts  what  it  needs  from  the  data¬ 
base  and  returns  wfat  it  concludes  to  the  database. 
Modules  may  then  interact  w  ith  the  results  of  other 
modules  by  examining  data  stored  in  appropriate 
places. 

Seismic  data  are  stored  in  a  manner  like 
that  used  at  the  Center  for  Seismic  Studies,  i.e.  as 
w  av  eform  files  and  a  set  of  associated  tables  w  hich 


describe  each  waveform,  but  with  enhancements. 
The  tables  used  in  YES  are  a  superset  of  those 
defined  in  the  CSS  version  3.0  database  specifica¬ 
tion  (Anderson,  cru/.,  1990).  Among  the  additions 
made  at  S-CUBED  are  tables  which  allow  the 
simultaneous  presence  of  several  sets  of  observed 
and  derived  information,  along  with  a  facility  to 
specify  which  of  the  several  sets  contains  the 
currently  preferred  information.  This  feature’s 
primary  purpose  is  to  impart  traceback  capability, 
so  that  an  investigator  can  determine  why  yield 
determinations  of  an  event  made  at  several  ses¬ 
sions  differ.  For  example,  as  more  stations  report 
their  locations  for  an  event,  the  origin  information 
associated  with  a  set  of  waveforms  ma>  change; 
previous  determinations  of  >ield  may  have  de¬ 
pended  upon  the  origin  as  it  w  as  know  n  at  the  time, 
and  recreation  of  those  previous  results  depends 
upon  recovery  of  the  information  about  origin 
known  at  the  time.  Station  and  event  magnitudes, 
station  corrections,  and  input  to  analysis  programs 
are  all  saved  in  a  way  that  permits  later  traceback. 
Waveforms  viewed  with  the  analyst’s  station  al¬ 
ways  use  the  currently  preferred  set  of  waveform 
tables. 

The  database  is  designed  to  facilitate  data 
flow  between  YES  analysis  modules.  In  general, 
the  procedure  used  b>  analysis  programs  is  to 
extract  waveform  and/or  tabular  data  from  the 
database,  measure  the  size  of  the  ev  ent,  and  append 
the  result  and  all  important  processing  parameters 
to  the  database.  The  analy  st  station  can  sav  e  arriv  al 
times  and  magnitudes  measured  on  individual 
waveforms  to  the  database.  Some  of  the  magni¬ 
tude  measurement  modules  retrieve  magnitudes 
derived  from  individual  waveforms  and  merge 
them  to  obtain  a  network  magnitude,  while  others 


examine  waveforms  diiectl),  without  user  inter¬ 
vention,  then  average  the  results  to  obtain  a  single 
magnitude  measure.  The  yield  estimation  proce¬ 
dure  recovers  the  results  of  the  magnitude  mea¬ 
surement  modules  from  the  database  and  produces 
a  w  eighted  y  ield  for  the  ev  ent.  This  estimated  y  ield 
is  then  compared  with  yields  from  other  events  in 
the  database  to  derive  statistics  describing  the 
likelihood  that  the  ev  ent  under  examination  w  as  in 
compliance  with  or  in  contravention  to  existing 
treaties.  All  these  procedures  share  their  results 
via  the  database. 

2 . 4  Database  Contents 


As  mentioned  above,  the  foundation  of 
YES  is  a  comprehensive  database  that  allows  a 
user  to  extract  parameters  from  explosion  wave¬ 
forms,  combine  these  parameters  to  obtain  an 
estimate  of  the  explosion  yield,  and  review  these 
results  in  the  context  of  the  env  ironment  in  w  hich 
the  explosion  took  place.  Such  a  database  neces¬ 
sarily  contains  a  large  v  ariety  of  both  seismic  and 
non-seismic  data. 

The  seismic  component  can  be  broken  dow  n 
into  four  major  classes  of  information,  these  being; 

•  event  and  origin  (epicentral)  information 

•  w  av  eform  and  w  av  eform  related  information 

•  measurements  obtained  from  the  waveforms 

•  static,  tabular  information,  such  as  trav  el-time 
tables,  earth  models,  station  locations  and  in¬ 
strumentation. 

Each  of  these  major  classes  by  itself  can  be 


further  subdivided  with  finer  distinctions.  For 
example,  event  and  origin  infoimation  is  obtained 
from  many  classified  and  unclassified  sources, 
each  of  w  hich  prov  ides  the  information  at  different 
lev  els  of  detail.  Unclassified  locations,  origin  times 
and  body -wave  magnitudes  have  been  obtained 
from  the  ISC  and  NEIS  for  all  explosions.  Other 
origin  times  and  locations,  computed  using  a  joint 
epicenter  inv  ersion  technique,  hav  e  been  extracted 
from  Marshall  ct  al.  (1984).  Classified  seismic 
locations  and  magnitudes  computed  using 
USAEDS  network  stations  also  have  been  incor¬ 
porated  into  the  database,  as  have  locations  deter¬ 
mined  by  analysis  of  classified  and  unclassified 
satellite  photographs.  Event  yields  based  on 
AFTAC  analysis  of  USAEDS  waveforms  are 
available  for  all  events,  and  additional  yields  are 
available  from  the  published  Soviet  literature 
(Bocharov  ctal.,  1989),  from  data  exchanges  with 
the  Sov  iet  government,  and  in  the  case  of  the  JVE, 
by  direct  measurement. 

Sources  of  waveform  data  are  just  as  di¬ 
verse.  Wav  eform  data  are  received  routinely  from 
AFTAC,  CSS,  IRIS,  and  other  sources,  and  can  be 
further  categorized  by  type,  such  as  short-period, 
long-period,  broadband,  single  station,  array,  and 
3-component,  or  by  network,  such  as  USAEDS, 
NORSAR,  GDSN,  CDSN,  and  IRIS.  As  of  this 
date  the  waveform  database  contains  more  than 
1 0000  discrete  w  av  eforms  and  occupies  storage  of 
ov  er  300  megaby  tes  (MB ).  An  extensi  v  e  effort  has 
been  made  to  insure  high  quality  in  the  waveform 
databa.se.  Thus,  every  waveform  has  been  visually 
inspected  by  a  trained  analy  st  to  determine  w  av  e- 
form  quality.  In  the  course  of  reviewing  the  data 
many  of  the  collected  waveforms  were  found  to  be 
unusable  due  to  w  indow  ing  problems,  contamina 


tion  by  other  events,  timing  and  calibration  errors, 
spikes,  data  dropouts  and  other  instrument  prob¬ 
lems.  Thus,  it  was  found  that  the  data  preview 
effort  was  critically  important  in  the  development 
of  a  database  which  would  fulfill  the  goals  of  this 
project. 

Measurements  consist  of  those  produced 
by  the  analysis  modules  within  YES  itself,  as  well 
as  measurements  received  fiom  outside  sources, 
such  as  AFTAC  single-station  magnitude  mea¬ 
surements.  This  variety  is  necessary  for  perform¬ 
ing  detailed  comparisons  between  >ields  derived 
from  YES  and  those  computed  b>  other  methods 
using  different  data  sets. 

Tabular  information  consists  of  a  complete 
station  location  file  containing  both  classified  and 
unclassified  station  locations,  standard  travel-time 
tables  for  the  important  seismic  phases  (Herrin  et 
cil,  1968),  surface  wave  path  corrections  for  use  in 
moment  tensor  inversion,  and  a  comprehensive 
instrument  response  database.  The  instrument 
database  is  a  compilation  of  information  collected 
from  the  USGS,  AFTAC,  IRIS,  and  others  which 
details  the  characteristics  of  the  recording  instru¬ 
ments  as  a  function  of  station,  channel  and  date. 
There  are  currently  750  different  response  histo¬ 


ries  in  this  section  of  the  database. 

In  addition  to  seismic  data,  the  database 
contains  much  v  isual  and  geologic  data.  There  are 
approximately  44  MB  of  SPOT  panchromatic  im¬ 
age  data  for  the  Shagan  River  test  site,  some  of 
which  has  30m  resolution  and  is  used  to  view  the 
entire  test  site  at  once,  and  most  of  which  has  10m 
resolution  for  use  in  detailed  c.xamination  of  the 
neighborhood  of  a  particular  event.  Geologic  im¬ 
ages,  digital  topography ,  and  data  used  in  creating 
cross-sections  encompass  12MB.  About  1  MB  of 
geologic  contour  data  for  display  on  top  of  other 
images  were  derived  from  the  digital  geologic 
data.  The  geologic  data  were  provided  by  Dr. 
William  Leith  of  the  USGS,  based  on  on-site 
investigations,  treaty  exchange,  interpretation  of 
satellite  photos,  and  research  in  the  Soviet  litera¬ 
ture.  Topographic  information  came  from  the 
Defense  Mapping  Agency. 

This  combination  of  seismic  and  non-seis- 
mic  data  provides  an  environment  whereby  all 
av  ailable  information  regarding  a  particular  explo¬ 
sion  under  investigation  is  available  to  the  system, 
in  a  manner  which  permits  the  analyst  to  effec¬ 
tively  integrate  it  into  the  best  possible  yield  esti¬ 
mate  for  the  explosion. 


3  Overview  of  System  Capabilities 


The  overall  design  philosophy  which  has 
been  followed  in  the  implementation  of  the  YES 
has  been  described  in  Section  2  and  in  a  preceding 
annual  report  (Murphy,  1990).  In  this  section, 
some  of  the  capabilities  and  functionality  of  the 
bj  stem  w  ill  be  graphically  illustrated  through  dis¬ 
plays  of  the  screens  encountered  b>  an  analyst  in  a 
typical  processing  session  for  a  selected  explosion. 
For  this  example,  the  unclassified  data  from  the 
Soviet  JVE  explosion  of  14  September  1988  will 
be  analyzed  and  used  to  illustrate  various  features 
of  the  system. 

As  has  been  noted  previously,  a  distin¬ 
guishing  characteiistic  of  the  YES  is  that  it  is 
completely  menu-driven  and  mouse-activated  and 
requires  no  keyboard  entry  by  the  user.  The  top 
level  menu  prov  iding  access  to  the  system  is  show  n 
in  Figure  1  where  it  can  be  seen  that  it  consists  of 
six  “buttons”  which  can  be  used  to  initiate  (SITE, 
EVENT,  PHASE,  FUNCTION)  or  terminate 
(QUIT)  action  within  the  system  or  to  view  online 
information  regarding  the  operating  characteris¬ 
tics  and  parameters  of  the  system  (HELP).  Select¬ 
ing  any  of  these  buttons  with  the  mouse  causes  a 
series  of  pulldown  menus  to  be  activated  as  illus¬ 
trated  in  Figure  2.  In  this  example,  the  SITE  button 


has  been  activated  to  select  the  Soviet  Shagan 
River  (Balapan)  test  site,  and  the  EVENT  button 
has  been  activated  to  select  the  JVE  explosion 
w  hich  w  as  detonated  at  that  test  site  on  14  Septem¬ 
ber  1988.  Once  the  test  site  and  event  have  been 
selected,  the  remaining  analy  st  interaction  w  ith  the 
system  is  initiated  through  the  PHASE  and  FUNC¬ 
TION  buttons.  The  PHASE  button  provides  the 
analyst  with  the  capability  to  choose  from  among 
the  six  different  seismic  phases  listed  in  Figure  2 
for  which  digital  waveform  data  are  currently 
available  on  the  system.  The  FUNCTION  button 
provides  access  to  the  seven  principal  computa¬ 
tional  and  analy  sis  modules  w  hich  permit  the  ana¬ 
lyst  to: 

•  view  the  seismic  data  w  ithin  the  context  of  the 
available  information  regarding  the  specific 
test  location  under  investigation  {Satellite 
Image,  World  Map) 

•  interact  w  ith  the  recorded  seismic  data  to  pro¬ 
cess  it  and  extract  the  v  arious  magnitude  mea¬ 
sures  of  intere.>t  {Analyst  Station,  Magnitude 
Measurement) 


Figure  1. 

Main  menu 
structure  for 
YES. 
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Figure  2. 

Illustration  of  pulldown 
access  to  individual  menu 
functions. 


•  formally  combine  the  seismic  measures  of 
source  size  to  obtain  an  optimum  measure  of 
explosion  yield  and  iiuantitative  measures  of 
the  uncertainty  in  that  estimate  ( Yield  Estima¬ 
tion) 

•  statistically  assess  the  results  with  respect  to 
any  existing  treaty  thresholds  or  other  yield 
levels  of  particular  interest  {Statistical  Sum¬ 
mary,  Spreadsheet). 

Hav  ing  specified  a  test  site  and  a  particular 
explosion  (in  this  case  the  JVE  event),  a  typical 
analysis  sequence  would  begin  with  the  selection 
of  the  Satellite  Image  option  from  the  FUNC¬ 


TION  menu,  which  brings  to  the  screen  the  SPOT 
satellite  image  display  of  the  Shagan  Rivet  test  site 
shown  in  Figure  3.  In  this  initial  display  of  the  test 
site  information  interface,  the  locations  of  prev  i- 
ous  explosions  at  this  site  are  shown  as  color- 
coded  square  overlays,  with  the  current  event  high¬ 
lighted  by  a  yellow  diamond.  In  this  case,  the 
different  colors  are  used  to  differentiate  those 
explosions  about  w  hich  the  Sov  iets  hav  e  published 
data  in  the  open  literature  (blue)  from  those  for 
which  only  seismic  information  is  available  (red). 
It  is  important  to  note  that  this  is  not  merely  a  static 
display,  but  that  in  fact  the  image  and  the  overlays 
to  it  are  formally  lied  to  an  extensive  online  data¬ 
base  of  supplementary  information.  Thus,  foi 


Figure  3.  SPOT  satellite  image  of  the  Slw.gau  River  test  site  with  superimposed 

locations  of  the  historical  explosions  (squares)  and  current  event 
(diamond). 


'i 


example,  in  thib  figure  the  operator  hab  pointed 
with  the  moube  to  one  of  the  blue  bquareb  (arrow ), 
which  hab  initiated  a  procebb  b>  which  available 
information  about  that  event  hab  been  extracted 
from  the  datababe  and  dibpla>ed  on  the  informa¬ 
tion  line  below  the  menu  buttoiib,  indicating  that 
thib  explobion  wab  detonated  at  Shagan  on  10 
December  1 972  at  the  bpecified  latitude  and  longi¬ 
tude  and  that  Bocharov  et  al.  ( 1 989 )  hav  e  reported 
the  depth  ab  478m  and  the  yield  ab  140  kt. 

For  purpobCb  of  dibplay ,  the  .salellile  image 
bhown  in  Figure  3  hab  been  comprcbbcd  to  an 


effective  rebolution  of  about  30m  bo  that  the  entire 
tcbt  bite  can  be  viewed  on  a  bingle  bcreen.  How¬ 
ever,  the  data  correbponding  to  the  full  10m  rcbo- 
lution  of  the  SPOT  image  can  albo  be  viewed  by 
activ  ating  the  appropriate  button  (Full ReS)  on  the 
right  hand  margin  of  thib  dibplay  andbimply  point¬ 
ing  with  the  moubc  to  any  location  on  the  image. 
For  example.  Figure  4  bhovvb  the  full  rcbolution 
bub-image  correbponding  to  the  location  of  the 
cratering  explobion  of  15  January  1965  which 
dammed  the  Shagan  Rivei,  producing  the  promi¬ 
nent  lake  in  the  bouthcabt  qiuidraiU  of  Figure  3. 
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Pull  resolution  SPOT 
satellite  image  of  the  region 
surrounding  the  Shagan 
River  cratering  explosion  of 
15  January  1965. 


Another  feature  provided  by  this  image  display  Figure  5  illustrates  this  feature  by  way  of  a  corn- 
module  is  the  capability  to  interactively  adjust  the  parison  of  the  nominal  display  (right)  w  ith  that 
contrast  and  brightness  using  the  slider  bars  lo-  resulting  from  interactively  reversing  the  contrast 
cated  at  the  bottom  of  the  right  hand  menu  margin,  and  uccrcasing  the  brightness  (left)  by  reposition- 


Illustration  of  interactive  modification 
of  brightness  and  contrast  in  SPOT 
satellite  image  display. 


ing  thebe  blidei  b  w  ith  the  moube.  It  can  be  been  that 
different  featureb  are  emphabized  in  thebe  two 
dibplajb,  which  fiicilitateb  anal>bt  identification  of 
both  geologic  and  man  made  featureb. 

In  addition  to  the  event  locatioiib,  a  number 
of  geologic  and  topographic  datababCb  have  been 
“regibtered"  to  thib  SPOT  image  and  arc  available 
for  dibpla>  from  the  OVERLAYS  and  IM.AGES 
menub.  For  example.  Figure  6  bhowb  the  overla> 
of  Leith’b  (1989)  color-coded,  burfacc  geologic 
map  of  the  area,  where  the  location  of  the  current 
event  (bquare)  ib  .shown  in  the  context  of  the 
varioub  expo.sed  geologic  unit.s  ...id  the  trace  of  the 
prominent  Chinrau  fault  which  inteibeclb  this  por¬ 


tion  of  the  tcbtbite.  The  Defenbe  Mapping  Agenc>  'b 
topographic  databa.se  for  the  area  ib  dibpla>ed  in 
coloi -coded  image  form  in  Figure  7  together  with 
overlays  of  the  correbponding  topography  con- 
touib  and  current  event  location  .sclecied  from  the 
OVERLAYS  menu.  Alternately,  the  variation  of 
.surface  topography ,  as  well  as  bubburface  geology 
to  a  depth  of  about  1  km,  along  any  .specified  line 
acro.sb  the  image  can  bo  accessed  using  the  cross- 
seetion  button  ( Crossseu;  v ..  .ho .  ight  hand  mar¬ 
gin  of  the  image  display .  This  feature  is  illusliated 
in  Figure  8  vvheic  the  right  hand  panel  shows  the 
line  .selected  by  the  .inalyst  I  y  pointing  with  the 
mouse  to  two  arbitrary  points  on  the  SPOT  image 


Figure  7. 


Color-coded 
representation 
of  DMA  topo¬ 
graphic  data 
for  the  Shagan 
River  test  site 
with  superim¬ 
posed  top 
graphic  con¬ 
tours  and 
current  event 
location. 
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and  the  left  hand  panel  shows  the  resulting  vertical 
section  through  that  line  which  was  automaticall> 
produced  by  the  system,  together  with  the  loca¬ 
tions  and  approximate  depths  of  penetration  of 
explosion  emplacement  holes  encountered  along 
that  line.  As  with  the  surface  topography,  the 
variation  in  the  depth  to  any  selected  geologic 
interface  can  also  be  exhibited  in  image  forr.i,  as 
illustrated  for  the  granite  surface  in  Figure  9  where, 
once  again,  the  corresponding  depth  contours  and 
location  of  the  current  event  have  been  ov  erlaid  for 
reference  purposes.  It  is  immediately  e  v  ident  from 
this  presentation  that,  unlike  the  current  JVEevent, 


most  previous  explosions  at  this  test  site  have  not 
been  detonated  in  granite  in  that  the  granite  surface 
throughout  much  of  this  region  lies  uniformly 
below  typical  explosion  emplacement  depths. 

Having  completed  an  initial  rev  iew  of  the 
available  information  regarding  the  explosion  test 
environment,  the  analyst  can  next  proceed  to  an 
examination  of  the  corresponding  recorded  seis¬ 
mic  data.  As  an  initial  step,  an  overview  of  the 
locations  of  ;  'ations  for  w  hich  data  from  the  se¬ 
lected  explosion  are  av  ailable  on  the  sy  stem  can  be 
obtained  by  selecting  the  World  Map  option  from 
the  FUNCTIONS  menu.  The  resulting  displays 


Figure  9.  Color-coded  repi  estniation  of  depth  to  the  top  of  the  granite  surface 

beneath  the  Shagjn  River  test  site  with  superimposed  depth  contours 
and  current  event  location. 
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Hjgure  1 6.  World  map  projections  (A<  JOO^)  showing  locations  of  stations  for  which 
digital  teleseismic  P  (left)  and  Lg  (right)  data  are  available  for  the  current 
event. 
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for  the  JVE  teleseismic  P  and  regional  Lg  phases 
are  shown  in  Figure  10,  plotted  on  azimuthal 
equidistant  projections  of  the  globe  (A  <  100°) 
centered  on  the  Shagan  River  test  site.  As  illus¬ 
trated  here,  a  menu  item  (STATIONS)  within  this 
module  provides  the  user  with  the  capability  to 
select  any  one  of  these  recording  stations  from  a 
list  and  subsequently  cause  a  straight  line  to  be 
drawn  on  the  screen  between  the  .source  and  the 
receiver.  In  this  transformation,  the  straight  line 
corresponds  to  the  great  circle  path  and  thus  pro 
vides  the  analyst  with  a  view  of  the  surface  projec¬ 
tion  of  the  propagation  path  followed  by  the  main 
energy  groups  between  these  two  points. 

The  seismic  data  themselves  can  be  ac¬ 
cessed  by  returning  to  the  main  menu  (Figure  2) 


and  sequentially  selecting  a  seismic  phase  from  the 
PHASE  menu  and  Analyst  Station  from  the 
FUNCTION  menu.  Thus,  for  example,  specifying 
the  phase  Lg  in  this  manner  automatically  initiates 
aprocess  by  which  the  Lg  recordings  from  the  JVE 
explosion  are  extracted  from  the  online  database  of 
digital  waveform  data  and  displayed  on  the  screen 
in  the  format  shown  in  Figure  11.  In  this  and 
subsequent  seismogram  displays,  the  stations  are 
ordered  by  increasing  epicentral  distance  and,  for 
the  Lg  and  surface  waves,  the  data  are  plotted  as  a 
function  of  group  velocity.  It  can  be  seen  in  this 
figure  that  the  traces  are  marked  by  a  vertical  line 
at  a  group  velocity  of  3.5  km/sec,  indicating  the 
nominal  expected  onset  time  of  the  Lg  pha.se.  This 
display  permits  the  analyst  to  quickly  assess  the 
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qualit>  of  the  data  and  identif>  an>  prominent  the  Magnitude  Measurement  option  from  the 
characteristics.  Thus,  for  example,  it  can  be  seen  FUNCTION  menu  a.^  illustrated  in  Figure  12. 
from  Figure  1 1  that  the  Lg  signals  are  relati\el>  Selecting  RMS  Lg  from  this  list  of  five  currentl> 
weak  at  the  IRIS  station  KIY  (IKIV)  and  at  axailabic  magnitude  measures  initiates  a  series  of 
NORSAR  (NO  I  AO),  consistent  with 


the  expected  effects  of  propagation 
across  the  Caspian  Sea  (KIV)  and 
along  the  complex,  far-regional  path 
to  Norway. 

Haviim  verified  that  the  avail- 
able  Lo  data  are  suitable  for  further 
processing,  the  analyst  can  proceed 
to  magnitude  estimation  by  selecting 


operations  b>  which  the  Lg  waveforms  recorded 
from  the  JVE  event  are  automatically  processed 
using  algorithms  described  by  Ringdal  (1983)  to 
obtain  single-station  and  network- averaged  mea¬ 
sures  of  the  Lg  magnitude,  MLg-  Since  this  pro¬ 
cessing  takes  some  time  lo  accomplish,  it  is  run  in 
the  background,  thereby  permitting  the  analyst  to 
mov  e  on  to  other  tasks  w  hile  it  is  being  completed. 
In  the  present  example,  the  analyst  has  selected 
Teleseismic  Pfrom  the  PHASE  menu  and  Ana¬ 
lyst  Station  from  the  FUNCTION  menu,  result¬ 
ing  in  the  waveform  display  shown  in  Figure  13. 
The  P  wave  seismograms  are  plotted  here  as  a 
function  of  reduced  time,  extending  from  10  sec¬ 
onds  before  to  15  seconds  after  the  signal  onset 


time  which  is  denoted  by  the  vertical  line  seg¬ 
ments.  Generally,  there  are  more  recordings  than 
can  be  uisplayed  on  a  single  screen  in  this  format 
and,  in  such  cases,  the  analyst  can  readily  page 
forward  and  backward  through  these  multiple 
screens  using  the  NEXT  and  PREVIOUS  menu 
buttons  at  the  top  of  the  analy  st  station  display . 

An  extensiv  e  set  of  signal  processing  and 
analysis  capabilities  has  been  assembled  within 
the  analyst  station  module,  as  indicated  by  the 
TOOLS  menu  display  in  Figure  13.  These  tools 
permit  the  analyst  to: 

•  assign  waveform  quality  {Set  Quality), 

•  redefine  arriv al  times  {Pick  Arrival ), 


Figure  1 3.  Analyst  station  display  of  selected  vertical-component 
teleseismic  P  wave  signals  for  the  current  event. 
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•  meabure  amplitudcb  and  periodb  (Expand,  mp, 
Ms), 

•  filter  the  data  (Demean/Detrend,  Filtei), 

•  normalize  to  a  common  inbtniment  rcbponbc 
(New  Instrument)  and 

®  compare  a  selected  w  a\  eform  w  ith  w  a\  eforms 
recorded  at  that  station  from  prc\  ions  explo¬ 
sions  (Compare,  Clone). 

Thus,  for  example,  it  can  be  seen  from 
Figure  1 3  that  the  P  \v  a\  e  signal  recorded  at  station 
MAJO  in  Japan  (bottom  trace)  is  obscured  b> 
microseismic  noise  lying  outside  the  signal  pass- 
band.  Selecting  Filter  from  the  TOOLS  menu 


causes  the  interacti\e  filter  interface  to  be  dis 
played  as  shown  in  Figure  14,  where  the  analyst 
has  used  the  mouse  to  adjust  the  slider  bars  to 
define  low  and  high  frequency  cutoffs  of  0.5  and 
5.0  Hz,  respectively,  and  to  specify  a  second  order 
roll-off  outside  th^.’.  band  by  setting  the  corre 
spending  pole  numbers  to  2.  Subsequent  selection 
of  any  trace  with  the  mouse  causes  this  filter  to  be 
automatically  applied  and  the  result  displayed,  as 
shown  here  for  station  MAJO. 

Application  of  the  Co/npare  function  from 
the  TOOLS  menu  provides  another  illustration  of 
the  powerful  capability  for  interactive  analysis 
which  is  available  in  the  analyst  station.  For 


Figure  1 4.  Example  of  the  specification  and  subsequent  application  of  a 

bandpass  filter  to  the  data  of  Figure  J3. 
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example,  by  sequentially  selecting  this  function 
and  the  station  KONO  waveform  from  Figure  14, 
a  search  through  the  database  is  automatically 
initiated  to  identify  all  other  available  P  wave¬ 
forms  recorded  at  that  station  from  pre  v  ious  Shagan 
River  explosions.  Generally,  not  all  of  the  explo¬ 
sions  identified  by  this  search  are  of  equal  impor 
tance  for  comparative  purposes  and  selection  of  a 
meaningful  subset  is  facilitated  by  examining  their 
locations  in  the  context  of  the  SPOT  image  test  site 
information  interface,  as  indicated  in  Figure  15. 
Using  this  interface,  the  analyst  can  interactively 
select  specific  events  with  the  mouse,  as  indicated 


here  by  the  diamond  overlays,  on  the  basis  of 
proximity  to  the  current  event  or  some  other  crite¬ 
rion,  and  display  the  corresponding  P  waveforms 
in  a  time  expanded  analyst  station  mode  such  as 
that  of  Figure  16.  Here  the  current  JVE  event 
recording  is  shown  at  the  top  of  the  figure  and  the 
waveforms  from  the  selected  comparison  events 
are  displayed  beneath  it  in  the  order  in  which  they 
were  selected.  Also  illustrated  in  this  figure  is  the 
use  of  the  Clone  feature  from  the  TOOLS  menu 
which  permits  the  operator  to  select  any  trace  on 
the  screen  with  the  mouse  (in  this  case  the  JVE 
recording),  create  a  color-coded  (red)  copy  and 


Figure  1 5.  SPOT  locations  of  Shagan  River  explosions  recorded  at  station 
KONO  in  Norway.  The  diamond  symbols  denote  those  events 
selected  for  comparative  analysis. 
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Figure  1 6.  Comparison  of  the  P  wave  signal  recorded  at  station  KONOfrom 
the  current  event  (top  and  red)  with  the  signals  recorded  at  that 
station  frc’n  the  selected  events  of  Figure  15. 


c'ag  it  to  any  other  location  on  the  screen  to 
conduct  detailed  waveform  comparisons.  This 
capability  permits  the  analy  st  to  quickly  assess  the 
consistency  of  the  current  observation  with  past 
experience  and,  in  this  case,  to  conclude  that  the 
JVE  recording  at  KONO  is  quite  consistent  with 
previous  observations  at  that  station  from  other 
explosions  located  in  that  area  of  the  test  site. 

Once  the  initial  w  av  eform  rev  iew  has  been 
completed,  the  analy  st  can  proceed  to  magnitude 
estimation  using  the  Expand  and  mp  functions 
fromtheTOOLS  menu.  By  .sequentially  selecting 
Expand  and  any  waveform  from  the  analyst  sta¬ 


tion  display ,  a  magnified  reproduction  of  that  w  av  e- 
form  is  produced  as  shown  in  Figure  17.  Then,  by 
employ  ing  the  malfunction,  the  mouse  can  be  used 
to  r- jsition  a  rectangle  so  as  to  define  the  peak-to- 
pcak  amplitude  and  half  period  of  the  selected 
cycle  of  motion  as  shown  in  this  figure.  After  this 
process  has  been  repeated  for  each  usable  trace, 
the  resulting  amplitudes  are  input  to  the  mb  estima¬ 
tion  module  using  the  Magnitude  Measurement 
menu  of  Figure  12  where  they  are  converted  to 
ground  motion,  corrected  foi  cpiccntral  distance 
and  station  effects  and  logarithmically  averaged  to 
obtain  a  network-averaged  mb  value  and  associ- 
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Figure  17. 

Illustration  of  the 
interactive  determi¬ 
nation  of  the  ampli¬ 
tude  and  period 
values  used  to  define 
single  station  mb 
values. 
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ated  uncertainty.  The  resulting  individual  station 
mb  values  for  the  current  JVE  explosion  are  dis¬ 
played  in  Figure  18,  together  with  the  estimated 
netw'ork-averaged  value  of  6.012.  This  display 
nodule  also  provides  the  capability  for  the  analyst 
to  interactively  eliminate  questionable  data  points 
and  then  to  recalculate  the  network  average  using 
the  ADD/DELETE  and 
RECALC  menu  but¬ 
tons  shown  at  the  top  of 
the  figure.  This  pro¬ 
cess  can  be  continued 
until  a  final  stable  mb 
estimate  is  obtained,  at 
which  time  it  is  written 
to  the  database  using 
the  designated  SAVE 
menu  button. 

Note  from  Fig¬ 
ure  18  that  at  this  point 
the  system  has  notified 


the  analy  t  that  the  RMS  Lg  magnitude  estimation 
has  been  completed,  as  indicated  by  the  icon  in  the 
upper  left  hand  corner  of  this  display.  Selecting 
this  icon  with  the  mouse  produces  the  display  of 
indi V  idual  station  and  netw  ork-av  eraged  MLg  ^  al- 
ues  shown  in  Figure  19.  It  can  be  seen  from  this 
figure  that  the  estimated  network-averaged  MLg 
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Figure  18.  Q 

Comparison  of 
individual 
station  and 
network- 
averaged  mb 
magnitudes 
determined  for 
the  current 
event. 


Figure  19. 

Comparison 
of  individual 
statioitand 
network- 
averaged 
Mhg  magni¬ 
tudes  cleter- 
minedfor 
the  current 
event. 


on  this  display  and  re¬ 
turning  to  the  SPOT 
image  test  site  informa¬ 
tion  interface  to  over¬ 
lay  the  quantity  mb  - 
MLg(NORS)  at  the  cur¬ 
rent  event  location  as  in 
Figure  20,  where  con¬ 
tours  of  this  quantity  de¬ 
rived  from  results  of 
previous  explosions 
have  also  been  overlaid 
for  reference  purposes. 
It  can  be  seen  from  this 


value  of  5.969  is  somewhat  lower  than  the  corre-  display  that  while  this  magnitude  difference  is 
spending  mb  value  of  6.012  of  Figure  18,  and  it  typically  negative  in  the  northeast  quadrant  of  the 
would  be  natural  for  the  analyst  to  question  whether  test  site  (by  as  much  as  0.15  magnitude  units),  it  is 
there  is  any  significance  to  this  difference.  This  generally  positive  in  the  vicinity  of  the  current  JVE 
issue  can  be  addressed  by  selecting  a  representa-  event  location,  with  an  average  value  close  to  the 
tivestation(e.g.,NORS)fromtheRESIDUALmenu  observed  value  of  0.043.  Thus,  this  capability 


Figure  20. 

Comparison  of 
the  observed 
value  of  mb  - 
MLg(NORSAR) 
for  the  current 
event  with 
contours  repre¬ 
senting  ob¬ 
served  variation 
of  that  param¬ 
eter  for  previ¬ 
ous  events. 
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permits  the  analyst  to  quickly  conclude  that  the 
current  result  is  entirely  consistent  with  previous 
experience  in  this  area  of  the  test  site. 

The  long-period  surface  wave  data  can  be 
accessed  in  the  same  manner  as  the  Lg  and  P  phase 
data  described  above,  resulting  in  analyst  station 
displays  such  as  that  shown  in  Figure  21  for  the 
vertical  component  Rayleigh  w  av  e  phase.  It  can  be 
seen  from  this  example  that,  although  these  traces 
are  aligned  according  to  a  common  group  velocity 
scale,  the  surface  wave  signals  are  quite  variable 
and  difficult  to  correlate  from  trace  to  trace.  This 
is  at  least  partially  due  to  the  fact  that  the  response 


characteristics  of  the  long-period  instruments  used 
to  record  these  data  vary  significantly  between 
stations.  Thus,  for  example,  the  station  MAJO 
trace  shown  at  the  bottom  of  this  figure  was  re¬ 
corded  through  a  relatively  broadband  system, 
with  the  result  that  the  long-period  signals  of 
principal  interest  are  obscured  by  higher  frequency 
arrivals.  In  order  to  provide  a  common  basis  for 
comparison,  the  TOOLS  menu  includes  a  function 
{New  Instrument )  which  permits  the  analyst  to 
transform  the  data  to  that  which  would  have  been 
observed  if  the  same  instrumentation  had  been 
employed  at  each  station.  This  feature  is  graphi- 


FigUfe  21 .  Analyst  station  display  of  selected  long-period  Rayleigh  wave 
signals  for  the  current  event. 
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cally  illubtrated  in  Figure  22  where  the  analyst  has 
selected  the  nominal  SRO  long-period  response 
and  converted  all  of  the  traces  in  this  display  to  that 
response  by  simply  designating  them  with  the 
mouse.  It  is  ev  ident  from  this  example  that  such  a 
capability  greatly  facilitates  quantitative  compari¬ 
sons  and  evaluation  of  data. 

Once  the  available  long-period  Rayleigh 
and  Love  wave  data  have  been  previewed,  the 
analyst  can  return  to  the  Magnitude  Measure¬ 
ment  menu  of  Figure  12  and  select  Moment 
Tensor,  thereby  initiating  a  process  by  which  the 
observ  ed  long-period  data  are  formally  inv  erted  to 


obtain  a  surface  wave  measure  of  seismic  magni¬ 
tude.  While  this  computationally  intensive  pro¬ 
cess  is  running  in  the  background,  the  analy  st  ne.xt 
selects  Spectral  P  from  this  same  menu  and 
proceeds  with  an  interactive  determination  of  the 
network- averaged  P  wave  spectrum  correspond¬ 
ing  to  the  current  event.  The  spectrum  estimated  in 
this  process  can  be  inverted  to  obtain  either  an 
equivalent  spectral  magnitude  or  a  direct,  model- 
based  yield  estimate  (Murphy,  1989),  and  Figure 
23  show  s  the  menu  structure  for  the  latter  option, 
where  it  is  indicated  that  the  Mueller/Murphy 
granite  (M/M  Granite)  source  model  has  been 


Figure  22. 


Example  of  the  specification  and  subsequent  application  of  the 
instrument  response  normalization  feature  to  the  data  of  Figure 


21. 


Figure  23.  Menu  options  for  the 
estimation  of  network- 
averaged  P  wave 
spectra. 


selected  for  the  current  JVE  event.  The  resulting 
automatic  fitto  the  attenuation-corrected,  observed 
spectrum  is  shown  in  Figure  24  where  ii  is  indi¬ 
cated  that  a  theoretical  spectrum  corresponding  to 
a  yield  ( W)  of  1 1 9  kt,  a  pP  -  P  delay  time  (to)  of  0. 82 
seconds  and  a  pP/P  amplitude  ratio  (A)  of  0.17 
provides  the  best  overall  fit  to  the  data.  Alter¬ 
nately,  the  analyst  can  elect  to  interactively  change 
these  automatically  determined  model  parameters 
using  the  plus  and  minus  buttons  located  in  the  box 
in  the  upper  right  corner  of  this  display,  to  explore 
solutions  corresponding  to  goodness  of  fit  criteria 
not  incorporated  in  the  automatic  algorithm. 

The  icon  appearing  in  the  upper  left  corner 
of  Figure  24  indicates  that  the  surface  wave  mo¬ 
ment  tensor  inversion  processing  has  now  been 


Figure  24.  Comparison  of  normalized  observed  and  best-fitting  theoreti¬ 

cal  network-averaged  P  wave  spectra  for  the  current  event. 


Figure  25.  Comparison  of  the  inferred  snrfo  e  wave  moment  tensor  solution 
for  the  current  event  with  the  co,  responding  path  normalized, 
observed  Rayleigh  (left)  and  Love  (right)  wave  amplitude  data. 
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completed,  and  selecting  it  with  the  mouse  pro¬ 
duces  the  graphical  solution  summar)  of  Figure  25 
where  the  model  fits  to  the  Ra>leigh  (leftj  and 
Love  (right)  wave  data  corresponding  to  a  tectonic 
release  F  factor  of  0. 1 3  are  displayed.  It  is  informa- 
tiv  e  to  examine  this  solution  in  the  context  of  other 
available  information  regarding  the  test  environ 
inent.  This  can  be  accomplished  by  returningonce 
again  to  the  SPOT  image  test  site  information 
interface  and  successively  overlaying  the  current 
and  historical  moment  tensoi  solutions,  as  well  as 
the  map  of  surface  geologic  features,  as  shown  in 
Figure  26.  In  this  display,  the  concentric  circles 
represent  the  moment  tensor  solutions  where,  foi 
the  historical  events,  the  ratio  of  the  diameters  of 
the  red-to-blue  circles  is  equal  to  the  inferred 
tectonic  release  F  factors  and  the  line  segments 


interior  to  the  circles  are  parallel  to  the  associated 
strike  of  the  tectonic  component.  The  solution  for 
the  current  (JVE)  event  is  differentiated  here  by  an 
alternate  color  scheme  in  which  yellow  and  light 
blue  circles  replace  red  and  blue,  respectively.  It 
can  be  seen  from  this  display  that  the  solution  for 
the  current  event  is  quite  consi.stent  with  those 
previously  determined  for  nearby  historical  events. 
Moreover,  it  is  ev  ident  that  the  inferred  strikes  of 
the  tectonic  releases  triggered  by  the.se  explosions 
are  parallel  oi  subparallel  to  that  of  the  mapped 
Chinrau  fault.  Thus,  by  using  the  le.st  site  informa¬ 
tion  interface,  the  analyst  is  able  to  quickly  verify 
that  the  surface  wave  moment  tensor  solution 
which  was  obtained  foi  the  curient  event  is  consis¬ 
tent  with  both  the  solutions  from  previous  nearby 
explosions  and  w  ith  the  regional  tectonic  env  iron- 


Figure  26.  Comparison  of  the  surface  wave  moment  tensor  solution  for  the 
current  event  (yellow  and  light  blue  concentric  circles)  with  those 
for  nearby  Shagaii  River  explosions  ( red  and  dark  blue  concentric 
circles)  and  with  the  surface  geologic  map  of  the  area. 


ment  in  which  the  test  was  conducted. 

Having  determined  the  various  seismic 
measures  of  source  si^e  (including  a  regional  Pn 
magnitude  not  discussed  above),  the  anal>st  can 
proceed  to  estimate  the  corresponding  explosion 
yield  and  the  associated  uncertainty  in  that  esti¬ 
mate  by  selecting  the  Yield  Estimation  feature 
from  the  FUNCTION  menu  of  Figure  2.  This  y  icld 
estimation  module  provides  an  estimate  of  unified 
yield  based  on  multiple  magnitude  measures  and 
associated  “extremal  confidence  limits”  corre¬ 
sponding  to  different  sets  of  constraints  specified 
by  the  analyst.  The  confidence  limits  are  obtained 


as  the  solution  of  a  nonlinear  programming  prob¬ 
lem  in  w  hich  estimates  of  y  ield  for  a  training  set  of 
explosions  are  minimised  and  maximized  over  the 
space  of  admissible  parameters  (Rodi,  1 989,  Rodi 
andMurphy,  1990>.  Input  to  this  model  consi.sts  of 
multiple  network-averaged  magnitudes  foi  the 
current  event  and  a  nominal  mt)/yield  relation 
ba.sed  on  some  combination  of  data  analysis  and 
expert  opinion,  which  foi  explosions  at  the  Shagan 
River  test  site  is  taken  to  be  (Murphy.  1990). 

mb  =  4.45  +  0.75  log  W 


The  current  module  can  pro\  ide  estimates 
corresponding  to  five  different  sets  of  constraints, 
parameterized  b>  the  uncertaint>  in  the  average 
V  alue  of  mb  at  100  kt  (mb(lOO)),  the  uncertaint>  in 
the  slope  (b)  of  the  mb/>ield  relation  and  the  upper 
bound  on  the  absolute  value  of  the  correlation 
coefficients  amongst  the  yield  estimation  errors 
for  the  individual  magnitudes  (max  P).  The  menu 
for  the  >  ield  estimation  module  is  snow  n  in  Figure 
27  together  with  sample  output  corresponding  to 
the  application  of  Models  1  and  3  to  the  magnitude 
data  foi  the  current  JVE  event.  In  these  displays, 
the  resulting  unified  yield  estimates  (W)  and  asso¬ 
ciated  upper  bound  uncertainty  factors  (F)  are 
shown,  together  with  the  yield  estimates  obtained 
from  the  individual  magnitude  values.  Note  that 
the  Pn  yield  value  is  shown  here  as  an  open  circle 
to  denote  the  fact  that  it  was  not  used  in  the  formal 
computation  of  W  and  F  due  to  the  fact  that  an 


adequate  training  set  of  previous  values  does  not 
currently  o  ist  for  that  magnitude  measure.  It  can 
be  seen  that  for  these  two  cases,  the  unified  yield 
estimate  is  about  1 2 1  kt  and  the  uncertainty  bounds 
range  from  a  factor  of  1 .96  for  the  case  (Model  1 ) 
in  which  mb(lOO)  is  specified  as  5.95  ±  0.05,  to  a 
factor  of  1.66  for  the  case  (Model  3j  in  which 
mb(  100)  is  given  as  exactly  5.95.  In  Figure  28,  the 
Model  1  solution  from  Figure  27  (top;  is  compared 
with  the  result  of  rerunning  Model  1  without  the 
surface  wave  moment  tensor  magnitude  (i.e.,  by 
using  the  slider  bars  at  the  bottom  of  the  menu  box 
to  assign  a  w  eight  of  zero  to  W(Mo; ).  It  can  be  seen 
that  this  modification  results  in  a  significant  reduc¬ 
tion  of  the  F  factor  from  1 .96  to  1 .55,  reflecting  the 
relatively  uncertain  yield  estimation  capability  of 
the  surface  w  av  e  moment  as  determined  from  the 
analysis  of  the  training  set  of  multiple  magnitude 
data. 


Figure  27. 

Ma/iu  and 
sample  output 
for  the  unified 
yield  estima¬ 
tion  module. 
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Figure  28. 

Comparison  of 
unified  yield  esti¬ 
mates  (W)  and 
associated  uncer¬ 
tainties  (F)  ob¬ 
tained  with  (top) 
and  without  (bot¬ 
tom)  the  surface 
wave  moment 
tensor  magnitude. 
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Gi\en  an  optimum  beismic  Cbtimate  of 
unified  yield  and  the  associated  uncertaint>,  the 
analyst  can  next  proceed  to  assess  whether  this 

Figure  29.  Statistical  assessment 
menu. 
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estimate  is  consistent  with  the  150  kt  threshold  of 
the  Threshold  Test  Ban  Treaty  (TTBT)  by  select¬ 
ing  the  Statistical  Summary  option  from  the 
FUNCTION  menu.  This  module  provides 
access  to  a  variety  of  statistical  compliance 
assessment  tests  developed  by  Mission  Re¬ 
search  Corporation  (Gray  et  al.  1990),  as 
indicated  by  the  menu  displayed  in  Figure 
29.  In  this  example,  the  analyst  has  re¬ 
quested  a  summary  of  the  results  of  three 
single  event  compliance  tests,  producing  the 
graphical  summary  shown  in  Figure  30,  where 
the  red  areas  under  the  distribution  curves 
provide  a  measure  of  the  probabilitj  of  a 
violation  of  the  1 50  kt  threshold.  The  formal 
statistical  results  are  summarized  in  the  boxes 
in  the  upper  left  hand  corners  of  lhe.se  fig¬ 
ures,  and  it  can  be  .seen  that  the  results  of  each 
ofthc.se  three  tests  indicate  that  the  seismic 
data  are  consistent  with  a  yield  of  le.s.s  than 
150  kt  at  a  specified  false  alarm  rale  of  less 
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Figure  30. 
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than  2.5  percent.  The  resolving  power  of  such  tests  initial  negotiation  of  the  TTBT. 

can  be  illustrated  b>  considering  the  outcome  Having  completed  the  formal  anal) sis,  the 

which  would  have  resulted  if  these  same  seismic  Spreadsheet  option  from  the  FUNCTION  menu 


data  had  been  observed  from  a  test  below 
the  water  table  at  NTS.  In  this  case,  the 
Shagan  mb/yield  relation  would  be  re¬ 
placed  by  (Murphy,  1981): 

mb  =  3.94  +  0.81  logW 


Figure  31 .  Compliance  test  result  (Test  1) 
for  the  'scenario  in  which  the  current  event 
data  were  obsen’ed  front  an  explosion  below 
the  water  table  at  NTS. 


resulting  in  the  revised  Test  1  outcome 
shown  in  Figure  31.  It  can  be  seen  that 
under  this  hypothesis,  the  yield  estimate 
would  be  363  ktand  it  would  be  concluded 
that  this  explosion  violated  the  1 50  kt  limit 
of  the  TTBT.  Thus,  this  simple  example 
graphically  illustrates  the  importance  of 
the  test  site  magnitude  bias  effect  which 
has  received  such  intense  study  since  the 
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can  be  used  to  pro\ ide  a  less  formal  environment 
for  evaluating  the  seismic  >ield  estimate  for  the 
current  event  in  the  context  of  the  results  obtained 
from  analjses  of  previous  explosions  at  that  test 
site,  as  well  as  any  available  independent  yield 
calibration  data,  as  illustrated  in  Figure  32.  In  this 
display ,  the  four  indiv  idual  seismic  >  ield  estimates 
and  corii-.'ponding  unified  >ields  are  listed  for 
forty  of  the  largest  Shagan  Riv  er  explosions,  to¬ 
gether  with  selected  calibration  fields 
(W(Bocharov)),  taken  here  for  the  purposes  of 
illustration  to  be  the  three  high  yield  values  pub¬ 
lished  by  Bocharov  etal.  (1988).  The  ratios  of  the 
unified  seismic  to  calibration  yields  are  listed  in 


the  last  column,  together  with  the  average  offset 
(Avg  Bias)  and  imcertaint>  (F;  inferred  from  these 
limited  calibration  data.  For  the  nominal  magni- 
tude/>  ield  relations  listed  at  the  top  of  this  displa> , 
the  average  offset  is  a  factor  of  1.07  with  an 
associated  F  factor  of  1.25.  Figure  33  shows  the 
result  of  interactive!)  adjusting  each  magnitude/ 
)ield  relation  in  such  a  manner  that  the  individual 
average  uffsets  are  minimized.  It  can  oe  seen  that 
these  modifications  eliminate  the  average  bias  in 
the  unified  seismic  >  ield  estimate  (i.e.,  Av  g  Bias  = 
1.00)  and  reduces  the  F  factor  to  1.20.  This 
example  illustrates  how  the  spreadsheet  module 
provides  a  capabilit)  for  the  anal)st  to  rapidl) 


Figure  32.  Spreadsheet  summary  comparison  of  seismic  yield  estimates 

for  the  current  event  (09/1 4/SS)  with  those  obtained  for  se¬ 
lected  previous  Shagan  River  explosions. 
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Figure  33.  Spreadsheet  summary  illiistraiing  the  results  of  interactively 
modifying  the  designed  magnitude/yield  relations. 
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explore  the  effect.^  of  allemale  hjpolhcxc^  on  ihc 
.seismic  j  ield  estimation  proi,v.AS. 

The  sample  anaUsis  .session  dc.scribcd 
abo\e  has  prov  ided  an  overv  iew  of  some  of  the 
capabilities  which  are  currentlv  available  within 
the  Yield  E.stimation  System.  As  wasev  idem  from 
the  menu  structurc.s  .shown  in  the  graphical  dis¬ 


plays.  there  arc.  in  addition,  many  other  features 
which  were  not  c.xcrci.scd  here  in  order  to  hold  the 
dc.scription  to  a  manageable  length.  However,  the 
complete  system  is  currently  operational  at  the 
DARPA  CSS  where  it  can  be  e.\crcihcd  over  its 
entire  range  of  functionality. 
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4  Summary  and  Future  Plans 


4.1  Summary 


In  this  report  we  have  presented  a  brief 
summary  of  the  current  status  of  the  ongoing 
research  investigations  directed  toward  the  devel¬ 
opment  of  a  comprehensive  new  seismic  yield 
estimation  system  for  underground  nuclear  explo¬ 
sions.  More  specifically  ,  a  preliminary  prototype 
version  of  this  system,  which  is  currently  opera¬ 
tional  at  the  DARPA  CSS,  has  been  described  in 
detail  and  its  functionality  has  been  graphically 
illustrated  through  a  sample  application  to  the 
seismic  data  recorded  from  a  selected  explosion. 

The  software  system  design  criteria  were 
reviewed  in  Section  2,  where  the  characteristics  of 
the  data,  analysis  tools,  database  relations  and 
graphical  user  interface  were  described  in  the 
context  of  their  integration  into  a  comprehensive 
system  for  seismic  yield  estimation  and  compli¬ 
ance  assessment.  This  discussion  included  an 
overview  of  the  conceptual  model  for  the  system 
and  provided  a  description  of  how  the  prototype 
version  has  been  implemented  in  a  Sun  color 
workstation  en  vironment  using  software  built  upon 
the  framework  of  the  X  Window  graphics  and 
Oracle  database  management  systems. 

This  was  followed  in  Section  3  by  a  dem¬ 
onstration  of  sy  stem  capabilities  in  which  a  com¬ 
plete  processing  session  was  graphically  illus¬ 
trated  using  data  recorded  from  the  Soviet  JVE 
explosion  of  14  September  1988.  In  this  presenta¬ 


tion,  reproductions  of  the  actual  workstation  screens 
encountered  by  an  analyst  in  such  a  session  were 
used  as  a  framework  for  describing  the  simple 
menu  structure  on  which  the  graphical  user  inter¬ 
face  has  been  developed.  Among  other  features, 
this  demonstration  illustrated  the  manner  in  which 
available  test  site  information  is  presented  to  the 
analyst  in  the  context  of  a  SPOT  satellite  image  of 
the  test  site,  in  a  format  which  permits  the  analyst 
to  interact  digitally  w  ith  the  image  in  a  w  orkstation 
environment  to  create  overlays  to  and  cross-sec¬ 
tions  through  the  image  to  display  event  locations, 
topography  ,  surface  and  subsurface  geologic  data 
and  a  variety  of  geophysical  parameters  of  poten¬ 
tial  interest  in  yield  estimation  analy  sis.  In  addi¬ 
tion,  this  sample  session  was  used  to  illustrate  the 
capabilities  provided  by  the  system  which  permit 
the  analyst  to  interact  with  the  recorded  seismic 
data  to  process  it  and  extract  the  v  arious  magnitude 
measures  of  interest,  to  formally  combine  these 
seismic  measures  of  source  size  to  obtain  an  opti¬ 
mum  measure  of  explosion  yield  and  quantitative 
measures  of  the  associated  uncertainty  and  to  sta¬ 
tistically  assess  the  consistency  of  this  seismic 
yield  estimate  with  any  existing  treaty  thresholds 
or  other  yield  levels  of  particular  interest. 

4.2  Future  Plans 


Now  that  a  preliminary  working  prototype 
has  been  successfully  implemented,  the  research 
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effort  is  beginning  to  focus  on  the  extension  of  the 
system  to  encompass  data  from  a  wider  range  of 
explosion  test  sites  and  stations  and  on  the  devel¬ 
opment  of  a  semi-automated  documentation  mod¬ 
ule.  More  specifically,  the  following  system  modi¬ 
fications  are  planned  for  the  next  phase  of  system 
development: 

•  The  system  will  be  extended  to  incorporate 
explosions  at  the  Soviet  Novaya  Zemlya  test  site. 
This  addition  will  include  a  comprehensive  seis¬ 
mic  database  consisting  of  all  available  digital 
seismic  data  recorded  from  these  explosions  at 
stations  of  the  US  AEDS,  GDSN,  NORS  AR,  CDSN 
and  IRIS  networks,  as  well  as  a  test  site  informa¬ 
tion  interface  based  on  a  SPOT  satellite  image  of 
Novaya  Zemlya. 

•  A  new  module  will  be  added  to  the  YES 
which  will  be  designed  to  permit  the  analyst  to 


validate  and  process  the  regional  seismic  data 
recorded  at  the  Designated  Seismic  Stations  (DSS  j 
provided  for  in  the  1990  Protocol  to  the  TTBT. 

•  The  statistical  y  ield  estimation  and  assess¬ 
ment  modules  in  the  current  version  of  the  YES 
will  be  modified  to  incoiporate  any  supplemental 
information  resulting  from  CORRTEX  yield  esti¬ 
mation  activity  or  on-site  inspections  of  the  type 
provided  for  in  the  1990  Protocol  to  the  TTBT. 

•  A  semi-automatic  event  report  generation 
module  will  be  added  to  the  YES  which  will 
provide  comprehensive  graphical  and  text  docu¬ 
mentation  of  the  yield  estimation  analyses  con¬ 
ducted  for  any  selected  explosion. 
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